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Abstract 

Background: The transmission of the malaria parasite Plasmodium falciparum from the human to the mosquito is 
mediated by dormant sexual precursor cells, the gametocytes, which become activated in the mosquito midgut. 
Because gametocytes are the only parasite stages able to establish an infection in the mosquito, they play a crucial role 
in spreading the tropical disease. The human-to-mosquito transmission triggers important molecular changes in the 
gametocytes, which initiate gametogenesis and prepare the parasite for life-cycle progression in the insect vector. 

Results: To better understand gene regulations during the initial phase of malaria parasite transmission, we focused on 
the transcriptome changes that occur within the first half hour of parasite development in the mosquito. Comparison 
of mRNA levels of P. falciparum gametocytes before and 30 min following activation using suppression subtractive 
hybridization (SSH) identified 126 genes, which changed in expression during gametogenesis. Among these, 17.5% had 
putative functions in signaling, 14.3% were assigned to cell cycle and gene expression, 8.7% were linked to the 
cytoskeleton or inner membrane complex, 7.9% were involved in proteostasis and 6.4% in metabolism, 12.7% were cell 
surface-associated proteins, 1 1.9% were assigned to other functions, and 20.6% represented genes of unknown 
function. For 40% of the identified genes there has as yet not been any protein evidence. 
For a subset of 27 genes, transcript changes during gametogenesis were studied in detail by real-time RT-PCR. Of 
these, 22 genes were expressed in gametocytes, and for 15 genes transcript expression in gametocytes was increased 
compared to asexual blood stage parasites. Transcript levels of seven genes were particularly high in activated 
gametocytes, pointing at functions downstream of gametocyte transmission to the mosquito. For selected genes, a 
regulated expression during gametogenesis was confirmed on the protein level, using quantitative confocal 
microscopy. 

Conclusions: The obtained transcriptome data demonstrate the regulations of gene expression immediately following 
malaria parasite transmission to the mosquito. Our findings support the identification of proteins important for sexual 
reproduction and further development of the mosquito midgut stages and provide insights into the genetic basis of 
the rapid adaption of Plasmodium to the insect vector. 
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Background 

Up to date, the tropical disease malaria is one of the 
most devastating infectious diseases in the world, caus- 
ing 216 million new cases and approximately 655.000 
casualties each year [1]. Disease treatment and control 
measures are undermined by the spread of drug resis- 
tance in malaria parasites, particularly in populations of 
Plasmodium falciparum, the agent responsible for mala- 
ria tropica (reviewed in [2]). 

The transmission of malaria parasites from the 
human to the mosquito is mediated by specialized 
sexual precursor cells, the intraerythrocytic gameto- 
cytes. Maturation of P. falciparum gametocytes from 
stage I to stage V takes approximately 10 days, and 
during this period the gametocytes maintain a stable 
cell cycle arrest (reviewed in [3,4]). The mature game- 
tocytes circulate in the humans blood stream, but 
remain dormant until they are taken up by a blood- 
feeding mosquito. 

When entering the mosquito midgut together with the 
blood meal, the gametocytes become activated from the 
dormant stage by external stimuli, i.e. a drop in tempera- 
ture and the contact with the mosquito-derived molecule 
xanthurenic acid (XA) (reviewed in [5,6]). Gametocyte 
activation leads to rounding up of the cell, followed by 
parasite egress from the enveloping erythrocyte, which in- 
volves the rupture of two membranes, the parasitophorous 
vacuole membrane (PVM) and the erythrocyte membrane 
[7] (reviewed [8]). During gametogenesis the microgame- 
tocyte replicates its genome three times in order to pro- 
duce eight motile microgametes. Following the fusion of 
micro- and macrogametes a zygote forms and develops 
into an infective ookinete within the following 24 hours. 
The motile ookinete possesses an apical complex which 
enables it to traverse the midgut epithelium before settling 
down and forming an oocyst between epithelium and 
basal lamina (reviewed in [5]). 

During gametocytogenesis, P. falciparum expresses a 
new set of genes important for sexual development 
[9-14]. Ingestion by the blood-feeding mosquito again 
triggers molecular changes in the sexual stage parasites, 
which help adjusting the gametocytes to the insect and 
which on the one hand initiate sexual reproduction and 
further development of the parasite in the vector, on the 
other hand prepare the emerging gametes for the hostile 
environment of the mosquito midgut. Noteworthy, the 
midgut stages have to persevere outside a host cell for 
more than one day. During this time period, the cells are 
highly vulnerable to the aggressive factors of the gut, 
which among others include bacteria as well as human 
immune cells, antibodies and complement proteins 
present in the blood meal, and this exposure results in 
an approximate 1000-fold loss of parasite abundance 
[15] (reviewed in [5,6]). 



Gametocyte maturation and gametogenesis are par- 
ticularly accompanied by the coordinated expression of 
numerous adhesive surface proteins, including the EGF 
domain-containing proteins Pfs25 and Pfs28, the 6-cys 
proteins Pfs230 and Pfs48/45, and the LCCL domain- 
containing PfCCp proteins. It is noteworthy that the 
majority of these proteins can be divided into two clas- 
ses: One class of the adhesion proteins, including Pfs230, 
Pfs48/45 and the six PfCCp proteins, is expressed within 
the parasitophorous vacuole (PV) of the developing ga- 
metocyte, and the majority of these proteins assemble to 
multimeric protein complexes [16,17] (reviewed in [18]). 
The adhesion proteins are subsequently present on the 
gamete surface, but expression of these proteins usually 
ceases during fertilization. The expression of the second 
class of surface proteins starts following parasite trans- 
mission to the mosquito, as was shown for Pfs25 and 
Pfs28, and expression often persists until the ookinete has 
formed (reviewed in [5]). One reason for this sudden 
onset of protein expression following gametocyte activation 
in the mosquito midgut is the translational repression of 
messenger RNA (mRNA) encoding for some of these pro- 
teins. This was inter alia shown in the rodent malaria para- 
site P. berghei for the repression of Pbs25 and Pbs28 by the 
RNA helicase DOZI (development of zygote inhibited) as 
part of a ribonucleoprotein complex [19]. DOZI and its 
subsequently identified interaction partner CITH (CAR-I 
and fly trailer hitch) are considered repressors of maternally 
supplied mRNA important for ookinete development [20]. 

While astonishing advances have been made in previous 
years in unveiling the molecular basis of parasite trans- 
mission from the mammalian host to the mosquito, the 
currently available data are fragmented and major key 
players for gametogenesis and adaptation have not yet 
been identified. We thus aimed at investigating the 
changes in transcript levels during gametocyte activation 
in order to gain in-depth knowledge on the molecular 
switch-over that takes place in the parasite during trans- 
mission. We selected the suppression subtractive hybri- 
dization (SSH) technique to experimentally screen for 
genes differentially expressed during gametocyte activa- 
tion, because this PCR-based method amplifies induced 
genes while simultaneously suppressing house-keeping 
genes, and is thus particularly suited to discover new or 
unexpected genes [21]. Our data provide insights into the 
regulated gene expression of malaria parasites during 
transmission to the mosquito and enable the identification 
of novel proteins with essential functions for the mosquito 
midgut stages, which might be suitable targets for trans- 
mission blocking interventions. 

Results and discussion 

To determine changes in the transcriptome of malaria 
parasites during the initial phase of transmission to the 
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mosquito, we conducted a SSH on mRNA isolated from 
mature non-activated gametocytes and from gameto- 
cytes at 30 min after in vitro activation, which have com- 
pleted gametogenesis [7]. A subtracted cDNA library 
was constructed and the resulting cDNA was assigned 
to the respective plasmodial genes. 

We identified a total number of 126 genes, for which 
expression levels changed in the gametocytes during 
activation. The majority of genes can be assigned to six 
major ontology groups (Figure 1, Table 1): 22 genes 
(17.5%) have putative functions in signaling, 18 genes 
(14.3%) are assigned to processes linked to cell cycle and 
gene expression, 11 genes (8.7%) can be linked to the 
cytoskeleton or the inner membrane complex (IMC), 10 
genes (7.9%) have putative functions in protein traffick- 
ing, stabilization and degradation (proteostasis), eight 
genes (6.4%) are linked to general metabolic functions, 
and 16 genes (12.7%) are proteins of the cell surface or 
PVM. Furthermore, 15 genes (11.9%) can be assigned to 
other functions, including protein synthesis and process- 
ing. A total of 26 genes (20.6%) encode for proteins with 
unknown function, eight of which have sequences for 
transmembrane domains and five of which encode for a 
signal peptide. 

Roughly 60% (76/126) of the SSH genes have been 
detected in the currently available P. falciparum game- 
tocyte proteomes, while for 40% there is as yet no 
protein evidence (Table 2) [9,12,13,22]. Almost identical 
numbers apply to the 113 orthologs of the rodent 
malaria parasite P. berghei; the 13 P. falciparum genes 
without P. berghei orthologs include amongst others the 
three etramp genes, pfsl6 and resa. 

The P. berghei orthologs of the identified P. falciparum 
genes were compared with a global gametocyte transcrip- 
tome analysis of P. berghei DOZI and CITH gene deletion 
mutants; destabilized mRNAs of these mutants are candi- 
dates for translational repression in non-activated gameto- 
cytes [19,20]. A total of 11 genes (8.7%) were identified as 
candidates regulated by DOZI and CITH (Table 3); eight 
of those lack any gametocyte protein evidence. One of the 




Figure 1 Ontology groups of P. falciparum genes with changes 
in transcript expression during gametocyte activation. The 

changes in transcript abundance were determined by SSH between 
the mRNA of mature non-activated gametocytes and gametocytes 
at 30 min p.a. in vitro. A total number of 126 genes were identified 
(the numbers of genes per group are indicated). 



destabilized genes was identified as a single peptide hit 
while the remaining two destabilized genes had 19 and 74 
peptide hits in total. A number of 16 additional genes 
were destabilized in the combined DOZI and CITH 
datasets indicating they could also be under translational 
control (data not shown). The data comparison indicates 
that the majority of genes identified via SSH, which 
change in transcript expression levels during gametocyte 
activation, are not translationally controlled in non- 
activated gametocytes, while many such genes have been 
identified to function in zygote to ookinete transformation 
[23]. In conclusion, gene regulation appears to be import- 
ant for the early events of gametocyte activation, i.e. for 
gametogenesis and fertilization, while translational repres- 
sion particularly plays a role for the expression of proteins 
important for ookinete development. 

We used real-time RT-PCR analyses to confirm the 
transcript changes in gametocytes during activation for a 
subset of the identified genes and to determine, if these 
genes become up- or down-regulated in their expression 
during gametocyte activation. Total RNA was isolated 
from immature (stages III and IV) and mature (stage V) 
gametocytes of P. falciparum strain NF54 and from 
activated gametocytes at 30 min post- activation (p.a.). 
Furthermore, we isolated RNA from mixed asexual 
blood stages of the gametocyte-less P. falciparum strain 
F12. Initially, the synthesized cDNA of each sample was 
tested for its stage-specificity by diagnostic RT-PCR, 
using primers for the asexual blood stage gene ama-1 
[24,25] and for the gametocyte-specific gene pfccp2 [26]. 
The tests confirmed that no pfccp2 transcript was 
present in the F12 sample, while no ama-1 expression 
was detected in the purified gametocyte samples of 
strain NF54 (Figure 2 A, left). The absence of ama-1 
signals in the gametocyte samples and the absence of 
pfccp2 in the F12 asexual blood stage sample further 
showed that these were devoid of any contamination by 
genomic DNA (gDNA). An additional test for gDNA con- 
tamination was performed by using primers specific for 
the gene hdacl (histone deacetylase 1; PF3D7_0925700). 
In all parasite samples, i.e. in samples of F12 asexual blood 
stages as well as in NF54 immature, mature and activated 
gametocytes, hdacl transcript was present (Figure 2A, 
right), as shown by diagnostic RT-PCR. In sample prepa- 
rations lacking reverse transcriptase, on the other hand, 
no hdacl -specific PCR bands were detected. 

We chose 27 genes from the SSH analysis for com- 
parison of transcript abundance, with representatives from 
all ontology groups. The gene pfsl6 was selected as an 
internal control, because it is known to be highly ex- 
pressed in gametocytes throughout development, while it 
is absent in the gametes [7,27,28]. Furthermore two sexual 
stage-specific genes, actin II and pfs25, served as external 
controls. Pfs25 is expressed in vesicular structures during 
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Table 1 Genes of the Plasmodium falciparum gametocyte transcriptome as identified by SSH (comparing mature 
non-activated gametocytes and gametocytes at 30 min p.a.) 



Plasmodb ID 


Protein 


Predicted features 


pred. 


Signaling 








PF3D7_ 


02 14600 


serine/threonine protein kinase 




204.3 


PF3D7_ 


J 356800 


serine/threonine protein kinase ARK3 




475.7 


PF3D7_ 


_0302100 


serine/threonine protein kinase PfCLK-4/SRPK1 


Phosphorylation of SR splicing factors 


157.4 


PF3D7_ 


J 230900 


atypical protein kinase, RIO family 




111.5 


PF3D7_ 


J 238900 


calcium/calmodulin-dependent protein kinase, PK2 




58.9 


PF3D7_ 


_0624000 


hexokinase, HK 




55.3 


PF3D7_ 


_05 15300 


phosphatidylinositol 3-kinase, PI3K 




255.9 


PF3D7_ 


_1 1 13100 


protein tyrosine phosphatase, PRL 




25.1 


PF3D7_ 


J 403900 


phosphatase 




35.2 


PF3D7_ 


J 464600 


phosphatase 


SP 


170.5 


PF3D7_ 


.0210600 


conserved Plasmodium protein 


SP, calcium-lipid-binding domain CaLB 


53.7 


PF3D7_ 


J 2381 00 


calcyclin binding protein 


Calcyclin-binding domain 


266.4 


PF3D7_ 


J 442200 


GTP binding protein 




118.7 


PF3D7_ 


J 358600 


zinc finger protein 


Zinc-Finger domain 


1 18.3 


PF3D7_ 


_0927200 


zinc finger protein 


Zinc-Finger domain 


189.7 


PF3D7_ 


J 004300 


zinc finger, C3HC4 type 




131.7 


PF3D7_ 


J 433400 


conserved Plasmodium membrane protein 


7 TMs, PHD C4HC3 zinc-finger-like motifs 


682.5 


PF3D7_ 


J 463900 


conserved Plasmodium membrane protein 


SP, 7 TMs, EF hand calcium-binding domain 


127.1 


PF3D7_ 


J 233600 


asparagine and aspartate rich protein 1, AARP1 


Zinc-finger domain, Ring/U-Box 


146.3 


PF3D7_ 


.0530200 


triose phosphate transporter iTPT 


SP, 7 TMs 


59.8 


PF3D7_ 


_0303900 


conserved Plasmodium protein 


SP, phosphatidyl-ethanolamine-binding 
protein PEBP domain 


23.2 


PF3D7_ 


_08 18200 


14-3-3 protein 


Regulative protein, binds signaling enzymes 
like kinases or TM receptors 


30.2 


Cell cycle and gene expression 






PF3D7_ 


J 445900 


DEAD/DEAH box ATP-dependent RNA helicase 




60.0 


PF3D7_ 


J 01 9000 


eukaryotic translation initiation factor subunit elF2A 




280.8 


PF3D7_ 


J 41 7200 


NOT family protein 


ncy U Id UUl l Ul Ll a\ I bv_l I |JL 


519 8 


PF3D7_ 


J 103800 


CCR4-NOT transcription complex subunit 1, NOT1 


4 TMs 


392.0 


PF3D7_ 


_0622900 


ApiAP2 protein 


Tr^n^rrintinn fpirtnr 


237.4 


PF3D7_ 


.0420300 


ApiAP2 protein 


1 1 al 1 bv_l l|JUUI 1 laLLUI 


400 2 


PF3D7_ 


_06 15400 


ribonuclease 




323.2 


PF3D7_ 


_08 11300 


CCR4-associated factor 1, CAF1 




200.2 


PF3D7_ 


J 034900 


methionine-tRNA ligase 




101.2 


PF3D7_ 


J 469700 


conserved Plasmodium protein 


Homology to the transcriptional mediator 
subunit Med6 


24.6 


PF3D7_ 


.0513600 


deoxyribodipyrimidine photolyase 




129.2 


PF3D7_ 


J 105500 


centrin-4 (CEN4) 




20.1 


PF3D7_ 


J 01 8500 


PHF5-like protein 


mRNA splicing factor RDS3 domain 


12.5 


PF3D7_ 


_0827800 


SET domain protein SET3 


Chromatin regulator, nuclear protein, 
in male gametocytes 


283.6 


PF3D7_ 


J 11 5200 


SET domain protein, SET7 


Chromatin regulator 


94.3 


PF3D7_ 


J 360900 


polyadenylate-binding protein 


RNA-binding domain 


44.8 
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Table 1 Genes of the Plasmodium falciparum gametocyte transcriptome as identified by SSH (comparing mature 
non-activated gametocytes and gametocytes at 30 min p.a.) (Continued) 



PF3D7J011800 QF122 antigen 

PF3D7J 347500 DNA/RNA-binding protein Alba 4, ALBA4 



KH domain type 1, RNA binding 



131.6 
42.2 



Cytoskeleton, IMC Motor complex 

PF3D7_0422300 alpha-tubulin II 

PF3D7J 246200 actin I 

PF3D7_0503400 actin-depolymerizing factor 1, ADF1 

PF3D7J 020200 kinesin 

PF3D7_0724900 kinesin-like protein 

PF3D7J 246400 myosin A tail domain interacting protein, MTIP 

PF3D7J251200 coronin 

PF3D7_0708000 cytoskeleton-associated protein 

PF3D7_05 15700 glideosome-associated protein 40, GAP40 

PF3D7_0918000 glideosome-associated protein 50, GAP50 

PF3D7J351700 Alveolin, ALV6 



Cytoskeletal element 
IMC component 



IMC component 
Actin-binding protein 
Cap-Gly domain 
10 TMs, IMC component 
SP, 2 TMs, IMC component 
IMC component 



49.7 

41.2 

13.7 

99.1 

222.6 

23.5 

69.0 

103.6 

51.8 

44.6 

151.2 



Proteostasis (Protein trafficking, stabilization and degradation) 



PF3D7_0807500 
PF3D7J 466300 
PF3D7_0815800 



proteasome subunit alpha type ( 
proteasome subunit RPN2 
vacuolar sorting protein, VPS9 



PF3D7_1 136400 conserved Plasmodium protein 

PF3D7_0215400 conserved Plasmodium protein 

PF3D7_081 8900 heat shock protein 70, HSP70-1 

PF3D7_091 7900 heat shock protein 70, HSP70-2 

PF3D7_0306200 activator of Hsp90 ATPase, AHA1 

PF3D7_0708400 heat shock protein 90, HSP90 

PF3D7_01 02200 ring-infected erythrocyte surface antigen, RESA 



Component oft he proteasome core particle 

Proteasome non-ATPase regulatory subunit 1 

Ortholog of Saccharomyces cerevisae 
vacuolar sorting protein 9 

TPR domain, mediates protein interactions 
and multi-protein complexes 

WD40-domain protein, assembly and 
stabilization of MPCs 



PHISTand DnaJ domains 



29.5 

132.8 

214.7 

107.1 

110.8 

73.9 
72.4 
41.6 
86.2 
126.5 



Metabolism, general 

mal_mito_2 cytochrome c oxidase subunit I, coxl 

PF3D7J 444800 fructose 1,6-biphosphate aldolase 

PF3D7_0806800 vacuolar proton translocating ATPase subunit A 

PF3D7_0922200 S-adenosylmethionine synthetase, SAMS 

PF3D7_0801900 conserved Plasmodium protein 

PF3D7_0403200 conserved Plasmodium protein 

PF3D7_051 1800 myo-inositol 1 -phosphate synthase 

PF3D7_01 04400 4-hydroxy-3-methylbut-2-enyl diphosphate reductase, 
LytB 



cell surface and PVM 

PF3D7J 035300 
PF3D7J 335300 

PF3D7J 228600 



glutamate-rich protein; GLURP 
reticulocyte binding protein 2, RH2b 

merozoite Surface Protein 9, MSP-9 



SP, 12 TMs, component of the respiratory 
chain in mitochondria 

Glycolysis pathway 

6 TMs 



FAD/NADP-binding domain 
Possibly acetyl-CoA-synthetase 



1 TM, protein export domain, non- 
mevalonate pathway 



SP, on parasite surface 

1 TM, on parasite surface, involved in 
erythrocyte binding 

SP, on merozoite surface, involved in 
erythrocyte binding 



53.0 

40.1 

123.0 

44.8 

552.0 

69.4 

69.1 

62.5 



141.1 
382.9 
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Table 1 Genes of the Plasmodium falciparum gametocyte transcriptome as identified by SSH (comparing mature 
non-activated gametocytes and gametocytes at 30 min p.a.) (Continued) 



PF3D7_ 


J 028700 


merozoite TRAP-like protein, MTRAP 


SP, 1 TM, belongs to TRAP family, in 
micronemes and on merozoite surface 


58.1 


PF3D7_ 


_0731500 


erythrocyte binding antigen 175, EBA-175 


SP, 1 TM, Duffy binding domain, on 
merozoite surface 


174.6 


PF3D7_ 


_01 02500 


erythrocyte binding antigen-181, EBA-181 


SP, 2 TMs, Duffy binding domain, on 
merozoite surface 


181.2 


PF3D7_ 


J 035700 


duffy binding-like merozoite surface protein, DBLMSP 


SP, Duffy binding domain, SPAM 


80.3 


PF3D7_ 


_0202500 


early transcribed membrane protein 2, ETRAMP2 


SP, 1 TM, integral PVM protein 


11.5 


PF3D7_ 


_0423700 


early transcribed membrane protein 4, ETRAMP4 


SP, 2TMs, integral PVM protein 


10.2 


PF3D7_ 


J 033200 


early transcribed membrane protein 10.2, ETRAMP 10.2 


SP, 2TMs, integral PVM protein 


38.9 


DE3P>7 

rrjU/_ 


_ I I UodUU 


LrVV-vVrL Tamny protein 


SP, adhesin, cysteine-rich CPW-WPC domain 


o/.o 


PF3D7_ 


_0406200 


Dfr 1 

rtS 1 0 


SP, 2 TMs, integral PVM protein of 
gametocytes 


1 ~1 
I 0./ 


PF3D7_ 


_0930300 


merozoite surface protein 1, MSP1 


SP, GPI anchor, EGF domain, on merozoite 
surface, part of MSP complex 


195.7 


PF3D7_ 


J 21 8800 


conserved Plasmodium protein, PSOP17 


SP, lectin domain, laminin G2 domain 


39.6 


PF3D7_ 


_0707300 


rhoptry-associated membrane antigen, RAMA 


GPI anchor, involved in erythrocyte binding 


103.6 


PF3D7_ 


J 021 900 


10b antigen 


Ankyrin repeat region, PHAX RNA-binding 
domain, interacts with surface proteins like 
PfEMPl'(Y2H) 


267.4 


Others 








PF3D7_ 


J 103100 


60S ribosomal protein P1 




13.0 


PF3D7_ 


J 027800 


60S ribosomal protein L3 




44.2 


PF3D7_ 


J 408600 


40S ribosomal protein S8e 




25.1 


PF3D7_ 


J 358800 


40S ribosomal protein S13 




17.3 


PF3D7_ 


.0206200 


metabolite/drug transporter 


1 1 TMs, Major facilitator superfamily 


62.7 


PF3D7_ 


J 207700 


blood stage antigen 41-3 precursor 


SP 


43.4 


PF3D7_ 


J 457000 


signal peptide peptidase, SPP 


8 TMs 


47.6 


PF3D7_ 


_0320700 


^innal npntiH^KP rnmnlpy Mihunif 1 SPC? 

-JlUI ICll |*JC|*J 11 ^1 ClOC Lul 1 l|*JICA omkJLJI IIL Z_, .Jl V_-Z_ 


2 TMs, SPC25 domain 


21 .0 


PF3D7_ 


_0207700 


cprinp rpnppt pnfinpn A ^FRA-4 
oCIII IC IC|JCdl dl lliyci 1 '-r, OCrvrA '-t 


^P n^npin fpmilv r\/c:fpinp nrnfpin^QP 
or, |jci|jciiii idiiuiy Lyjicnic pi uicn idoc 


1 08 7 


PF3D7_ 


J 033800 


nl a empnei n \/ll 
LJIdol 1 Icpoll 1 VII 


QP 1 TA/l acnarf\/l pirnfpacp 
Or, I I /VI, dj|Jdl lyi piULcdoc 


52 3 


PF3D7_ 


_0930000 


nrnrnll^npn l\/c:inp S-Hinvx/npn^p 
ui ului lay ci i iyoii ic .J uiUAyyci laoc 


2 TMs 


66.8 


PF3D7_ 


.0817600 


rnn^prx/prl Pl^^mnHii im nrnfpin 
^.ui i jci vcu r la ji i iuuiui i i piui_cni 


MAC/PF Hnmain 


1 1 2.2 


PF3D7_ 


_0403800 


fnn^pr\/pH PI^KmnHii im nrntpin 

Lul IjCI VCU 1 IQOIIIUUIUIII Ul ULCM 1 


Alnhs/hpfs-hvHrnls^p 

/\ikji iq/ ucia i ly ui uia jc 


83.4 


PF3D7_ 


_0709900 


rnnc:pr\/pH Pl^mnHii im mpmhnnp nrnfpin 

L.UI IjCI VCU r la jl 1 IUUIUI 1 1 1 1 ICI 1 1 Ul CM IC UIULCIII 


A TMq plnh^/hpfp-hvHrnl^QP 
i l ivio, dl|Jlld/UCld liyuiUldoC 


2740 


PF3D7_ 


_03 11600 


rnnc:pr\/pH Pl^mnHii im nrnfpin i inknnwn fi inrfinn 
L.UI i oci vcu r ia ji i iuu i ui 1 1 uiuiciii, ui im iuvvi i i ui iluui i 


^P 1 Tl\/I ri hnnhnri n-l Hnm^in ni ifpfi\/p 
jr, i i ivi, iiuuuiiuiiii i uui i id 1 1 i,uu Ld li vc 

dolichyl-diphospho-oligosaccharide-protein 


84 0 


Unknown 








PF3D7_ 


_0904200 


conserved Plasmodium protein 


SP, 1 TM 


34.0 


PF3D7_ 


J 024800 


conserved Plasmodium protein 


SP, 2 TMs 


171.1 


PF3D7_ 


J 239400 


conserved Plasmodium protein 


SP 


22.0 


PF3D7_ 


_0730800.1 


Plasmodium exported protein 


SP 


32.1 


PF3D7_ 


J 3 16700 


conserved Plasmodium protein 


SP 


73.4 


PF3D7_ 


_0601900 


conserved Plasmodium protein 


1 TM 


14.7 


PF3D7_ 


_0704100 


conserved Plasmodium membrane protein 


6 TMs 


425.3 


PF3D7_ 


J 023000 


conserved Plasmodium membrane protein 


5 TMs 


48.9 


PF3D7_ 


J 439600 


conserved Plasmodium protein 


1 TM 


22.6 
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Table 1 Genes of the Plasmodium falciparum gametocyte transcriptome as identified by SSH (comparing mature 
non-activated gametocytes and gametocytes at 30 min p.a.) (Continued) 



PRD7 1071700 


rnrKPivprl Pl^^mnHii im mpmhranp nrnfpin 
^.^jiiDCrivcru riaoiiiwuiuiii iiicriiikjiaiicr lji^jlctmi 


4 TMs 


830.2 




rnn^pn/pH Pl^mnHii im nrnfpin 


i TUc intprarts with PRD7 1071700 (Y7H) 

z. 1 1 vio, iii lci a^.Lo vviiii ri ulj i i wz. i / uu v i z.i \j 


8224 


PRD7 1076600 


rnn^prx/prl Pl^^mnHii im nrnfpin 
100 vcu r iaoi i iuuiui i i pi^LOi i 


Infpr3rf<; wifh RN A-^^nri^fprl nrnfpin<; (Y?Fn 

IIILCia^-LO Willi IMN/\ ajjULIULCU piUlClllj \ i 1/ 


196.5 


r 1 JLy/ 1 JvJZ-VjUU 


rnn^p^/pH Pl^mnHii im nrnfpin 
lui i jci vcu r la oi i iuu i u i i i piwicni 


Infpr^rfc; \A/ifh c;nlirinn farfnr ^iihiinif 
ii iicia^io vvili i o|-jii^.ii ly ia^.Lkji ju oukjui ml 

2 (Y2H) 


32.1 


PF3D7_0202400 


conserved Plasmodium protein 


Pfg27 domain 


141.9 


PF3D7_1 124200 


conserved Plasmodium protein 


UAS domain 


49.3 


PF3D7_1 452400 


conserved Plasmodium protein 




107.5 


PF3D7_1 225600 


conserved Plasmodium protein 




93.1 


PF3D7J451200 


conserved Plasmodium protein 




176.2 


PF3D7_1 438800 


conserved Plasmodium protein 




84.0 


PF3D7_1 114600 


conserved Plasmodium protein 




39.6 


PF3D7_1026100 


conserved Plasmodium protein 




23.9 


PF3D7J230100 


conserved Plasmodium protein 




61.5 


PF3D7_1321000 


conserved Plasmodium protein 




80.6 


PF3D7_0411000 


conserved Plasmodium protein 




179.7 


PF3D7_1 126700 


conserved Plasmodium protein 




110.1 


PF3D7_0508900 


conserved Plasmodium protein 




370.2 



IMC, inner membrane complex; PVM, parasitophorous vacuole membrane; SP, signal peptide; TM, transmembrane; Y2H, yeast-two-hybrid. 



gametocyte maturation and relocates to the surface of 
macrogametes following activation. Pfs25 is subsequently 
present on the parasite surface until the ookinete stage 
[16,29-32]. Actin II is a sexual stage-specific actin isomer, 
and in P. berghei actin II was reported to play a role 
during microgametogenesis [33,34]. In addition, the gene 
su a5, encoding for subunit (SU) aS of the a-ring of the 
proteasome core particle [35], was included in the 
investigations. 

Transcript expression levels were measured via real- 
time RT-PCR and calculated by the 2" ACt method [36,37] 
in which the threshold cycle number (Ct), was normalized 
to the Ct of the endogenous control gene encoding for 
P. falciparum seryl tRNA synthetase (PF3D7_0717700) as 
reference gene [38,39] . Transcript levels with 2" ACt values 
below 0.5 were considered negligible. Real-time RT-PCR 
revealed transcription in mature gametocytes for 22 out of 
the 27 genes. Out of these, 15 genes showed increased 
transcript expression in gametocytes compared to asexual 
blood stage parasites (Figure 2B). 

Seven genes were identified, i.e. PF3D7_0827800 (setS), 
Pf3D7_1246200 (actin 7), PF3D7_1218800 ipsoplT), 
PF3D7_1239400, PF3D7_0704100, PF3D7_0417400, and 
PF3D7_1321000, for which transcript levels were in- 
creased in activated gametocytes as compared to asexual 
blood stages, immature gametocytes and non-activated 
gametocytes (Figure 2B), indicating that these genes may 
play an important role downstream of gametocyte acti- 
vation in the mosquito midgut. SET3 was previously 



described to accumulate in male gametocytes, where it 
contributes to a prompt entry and execution of S/M 
phases during microgametogenesis [40]. SET domains are 
assigned to chromatin dynamics and are often found in 
histone methyltransferases, thus they play a role in the 
epigenetic control of gene regulation. The P. falciparum 
genome encodes for at least nine SET-domain-containing 
proteins which exhibit five different types of substrate 
specificities [41]. 

Actin I was described as part of the plasmodial motor 
complex (reviewed in [42]) and was recently also reported 
to be present in gametocytes [33,43]. Furthermore, 
PSOP17 (putative secreted ookinete protein 17) was previ- 
ously reported to be expressed in mature gametocytes of 
P. falciparum [9,12] and in ookinetes of P. berghei [44]. In 
addition, transcript expression of the two control genes, 
PF3D7_1031000 (pfs25) and PF3D7_1412500 (actin II) 
was strong in gametocytes compared to asexual blood 
stage parasites and increased during gametocyte activation 
(Figure 2B). 

Eight genes had an increased transcript expression in 
gametocytes compared to asexual blood stage parasites, 
and transcript levels either remained constant or decreased 
following activation. These include PF3D7_0302100 (clk4\ 
PF3D7_1115200 (setT), PF3D7_0422300 (a-tubulin II), 
PF3D7_1 103500 (encoding for a CPW-WPC protein), 
PF3D7_406200 (pjsl6), PF3D7_1457000 (spp), PF3D7_0817600 
(encoding for a MAC/PF domain, here termed PPLP6) 
and PF3D7_1225600 (Figure 2B). The mRNA splicing 



Table 2 Comparison of genes identified by SSH with P. falciparum and P. berghei gametocyte proteome data 



Pf P. falciparum GC P. berghei male/female GC Total ^ qj 



PFJD 


Gene Ontology 


MW 


Troph 


1 & II 


IV & V 


V 


I.II.IV.V 


FG1 


FG 3 


FG (1+3) 


MG1 


MG3 MG(1+3) 


FG +MG 


FG.MG.PF 


PBJD 


< Si 
1 


PF3D7_02 10600 


signaling 


53.7 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


PBANKA_030750 


tr do 
o" ^ 


PF3D7_02 14600 


signaling 


204.3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


PBANKA_031 140 


C Genomics 2013 
medcentral.com/ 


PF3D7_0302100 


signaling 


157.4 


2 


3 


1 


7 


1 1 


0 


0 


0 


0 


0 0 


0 


1 1 


PBANKA_0401 10 


PF3D7_0303900 


signaling 


23.2 


0 


4 


4 


1 1 


19 


0 


0 


0 


0 


0 0 


0 


19 


PBANKAJD40250 


PF3D7_05 15300 


signaling 


255.9 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


PBANKA_1 1 1490 


PF3D7_0530200 


signaling 


59.8 


0 


0 


0 


3 


3 


0 


0 


0 


0 


0 0 


0 


3 


P B AN KA_ 124460 


^ -» 


PF3D7_0624000 


signaling 


55.3 


17 


20 


2 


24 


46 


6 


13 


19 


1 1 


10 21 


40 


86 


PBANKA_1 12290 


:256 
1-21 


PF3D7_0818200 


signaling 


30.2 


18 


21 


14 


25 


60 


1 1 


20 


31 


9 


1 1 20 


51 


1 1 1 


PBANKA_071260 


64/1 


PF3D7_0927200 


signaling 


189.7 


0 


13 


2 


5 


20 


0 


0 


0 


0 


0 0 


0 


20 


PBANKA_082800 


4/256 


PF3D7_ 1004300 


signaling 


131.7 


3 


5 


1 


0 


6 


0 


0 


0 


0 


0 0 


0 


6 


PBANKA_1 20260 


PF3D7_1 1 13100 


signaling 


25.1 


1 


2 


1 


3 


6 


0 


0 


0 


0 


0 0 


0 


6 


PBANKA_093450 




PF3D7_1 230900 


signaling 


1 1 1.5 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


P B AN KA_ 144560 




PF3D7_1 233600 


signaling 


146.3 


0 


4 


0 


8 


12 


0 


0 


0 


0 


0 0 


0 


12 


PBANKA_1 44820 




PF3D7_1 2381 00 


signaling 


266.4 


7 


9 


3 


9 


21 


2 


6 


8 


8 


7 15 


23 


44 


PBANKA_1 45260 




PF3D7_1 238900 


signaling 


58.9 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


P B AN KA_ 145340 




PF3D7J 356800 


signaling 


475.7 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


PBANKA_1 13310 




PF3D7_1 358600 


signaling 


1 18.3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


PBANKA_1 13490 




PF3D7_1 403900 


signaling 


35.2 


6 


9 


5 


12 


26 


X 


X 


X 


X 


X X 


X 


26 


X 




PF3D7_1 433400 


signaling 


682.5 


0 


0 


2 


0 


2 


0 


0 


0 


0 


0 0 


0 


2 


PBANKA_101 150 




PF3D7_1 442200 


signaling 


1 18.7 


0 


0 


1 


0 


1 


0 


0 


0 


0 


0 0 


0 


1 


PBANKA_1 30610 




PF3D7_1 463900 


signaling 


127.1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


PBANKA_1 32710 




PF3D7_ 1464600 


signaling 


170.5 


23 


31 


0 


3 


34 


0 


0 


0 


0 


0 0 


0 


34 


PBANKA_1 32800 




PF3D7_0420300 


cell cycle & gene expr. 


400.2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


PBANKA_052170 




PF3D7_05 13600 


cell cycle & gene expr. 


129.2 


10 


29 


0 


10 


39 


0 


0 


0 


0 


0 0 


0 


39 


PBANKA_1 1 1330 




PF3D7_0615400 


cell cycle & gene expr. 


323.2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


PBANKA_1 23010 




PF3D7_0622900 


cell cycle & gene expr. 


237.4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


PBANKA_1 12180 




PF3D7_081 1300 


cell cycle & gene expr. 


200.2 


0 


2 


1 


2 


5 


0 


0 


0 


0 


0 0 


0 


5 


PBANKAJ 42620 




PF3D7_0827800 


cell cycle & gene expr. 


283.6 


0 


0 


1 


0 


1 


0 


0 


0 


0 


0 0 


0 


1 


PBANKA_070290 




PF3D7J011800 


cell cycle & gene expr. 


131.6 


41 


22 


6 


17 


45 


4 


11 


15 


14 


0 14 


29 


74 


PBANKAJ 21 020 


-n 


PF3D7_1018500 


cell cycle & gene expr. 


12.5 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


PBANKA_050270 


'age 


PF3D7_1019000 


cell cycle & gene expr. 


280.8 


7 


0 


0 


0 


0 


0 


0 


0 


0 


0 0 


0 


0 


PBANKA_050320 


00 
O 


PF3D7_1 034900 


cell cycle & gene expr. 


101.2 


21 


31 


0 


49 


80 


6 


0 


6 


3 


0 3 


9 


89 


PBANKA_051870 


— h 

NJ 



Table 2 Comparison of genes identified by SSH with P. falciparum and P. berghei gametocyte proteome data (Continued) 



v. $ 



PF3D7J 103800 


cell cycle & gene expr. 


392.0 


7 


12 


0 


23 


35 


0 


0 


0 


0 


0 


0 


0 


35 


PBANKA_094310 


% % 


PF3D7J 105500 


cell cycle & gene expr. 


20.1 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_094140 


i 


PF3D7J115200 


cell cycle & gene expr. 


94.3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_093250 


BMC 
bion 


PF3D7J 347500 


cell cycle & gene expr. 


42.2 


18 


15 


30 


22 


67 


11 


0 


11 


7 


10 


17 


28 


95 


PBANKAJ 36030 


CD Q 
Q_ fO 
n s 


PF3D7J 360900 


cell cycle & gene expr. 


44.8 


8 


11 


11 


12 


34 


X 


X 


X 


X 


X 


X 


X 


34 


X 


comics 2( 
entral.cc 


PF3D7J417200 


cell cycle & gene expr. 


519.8 


5 


6 


0 


6 


12 


0 


0 


0 


0 


0 


0 


0 


12 


PBANKAJ 02550 


PF3D7J 445900 


cell cycle & gene expr. 


60.0 


16 


25 


3 


21 


49 


0 


0 


0 


1 


0 


1 


1 


50 


PBANKAJ 30970 


)13, 


PF3D7J 469700 


cell cycle & gene expr. 


24.6 


0 


2 


0 


3 


5 


0 


0 


0 


0 


0 


0 


0 


5 


PBANKAJ 33290 


14:2 

471 - 


PF3D7_0422300 


IMC, cytoskeleton 


49.7 


11 


25 


13 


30 


68 


0 


0 


0 


0 


0 


0 


0 


68 


PBANKA_052270 


NJ <-n 
— » Q\ 


PF3D7_0503400 


IMC, cytoskeleton 


13.7 


6 


5 


1 


8 


14 


0 


0 


0 


0 


1 


1 


1 


15 


PBANKAJ 10310 


? 


PF3D7_05 15700 


IMC, cytoskeleton 


51.8 


3 


0 


0 


7 


7 


1 


0 


1 


0 


0 


0 


1 


8 


PBANKAJ 11 530 


NJ 
Ln 

ON 


PF3D7_0708000 


IMC, cytoskeleton 


103.6 


0 


0 


2 


12 


14 


0 


0 


0 


1 


0 


1 


1 


15 


PBANKA_080520 


PF3D7_0724900 


IMC, cytoskeleton 


222.6 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_062240 




PF3D7_0918000 


IMC, cytoskeleton 


44.6 


14 


12 


1 


23 


36 


0 


0 


0 


0 


0 


0 


0 


36 


PBANKA_081900 




PF3D7J 020200 


IMC, cytoskeleton 


99.1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_050430 




PF3D7_1 246200 


IMC, cytoskeleton 


41.2 


23 


28 


13 


29 


70 


12 


19 


31 


6 


6 


12 


43 


113 


PBANKAJ 45930 




PF3D7_1 246400 


IMC, cytoskeleton 


23.5 


1 


0 


0 


3 


3 


0 


0 


0 


0 


0 


0 


0 


3 


PBANKAJ 45950 




PF3D7_1251200 


IMC, cytoskeleton 


69.0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 46410 




PF3D7J351700 


IMC, cytoskeleton 


151.2 


0 


0 


1 


0 


1 


0 


0 


0 


0 


0 


0 


0 


1 


PBANKAJ 36440 




PF3D7_01 02200 


proteostasis 


126.5 


33 


3 


1 


0 


4 


X 


X 


X 


X 


X 


X 


X 


4 


X 




PF3D7_0215400 


proteostasis 


110.8 


0 


0 


1 


0 


1 


0 


0 


0 


0 


0 


0 


0 


1 


PBANKA_031210 




PF3D7_0306200 


proteostasis 


41.6 


4 


7 


0 


8 


15 


0 


0 


0 


0 


0 


0 


0 


15 


PBANKA_040460 




PF3D7_0708400 


proteostasis 


86.2 


54 


60 


17 


82 


159 


7 


0 


7 


8 


10 


18 


25 


184 


PBANKA_080570 




PF3D7_0807500 


proteostasis 


29.5 


11 


11 


2 


20 


33 


3 


8 


11 


4 


2 


6 


17 


50 


PBANKAJ 22310 




PF3D7_0815800 


proteostasis 


214.7 


0 


3 


0 


0 


3 


0 


0 


0 


0 


0 


0 


0 


3 


PBANKA_071500 




PF3D7_0818900 


proteostasis 


73.9 


45 


46 


33 


57 


136 


34 


40 


74 


31 


28 


59 


133 


269 


PBANKA_071190 




PF3D7_0917900 


proteostasis 


72.4 


45 


40 


42 


51 


133 


24 


46 


70 


16 


16 


32 


102 


235 


PBANKA_081890 




PF3D7_1 136400 


proteostasis 


107.1 


0 


3 


0 


11 


14 


1 


0 


1 


0 


0 


0 


1 


15 


PBANKA_091220 




PF3D7J 466300 


proteostasis 


132.8 


23 


26 


3 


34 


63 


4 


9 


13 


4 


0 


4 


17 


80 


PBANKAJ 32970 




mal_mito_2 


metabolism general 


53.0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_MIT0002 




PF3D7_01 04400 


metabolism general 


62.5 


0 


1 


1 


3 


5 


0 


0 


0 


0 


0 


0 


0 


5 


PBANKA_020870 


q; 


PF3D7_0403200 


metabolism general 


69.4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 00080 


ge 9 of 


PF3D7_05 11800 


metabolism general 


69.1 


36 


29 


2 


40 


71 


2 


0 


2 


0 


0 


0 


2 


73 


PBANKAJ 11 140 



NJ 



Table 2 Comparison of genes identified by SSH with P. falciparum and P. berghei gametocyte proteome data (Continued) 



v. $ 



PF3D7_0801900 


metabolism general 


552.0 


0 


0 


1 


0 


1 


0 


0 


0 


0 


0 


0 


0 


1 


PBANKAJ 22830 


% % 


PF3D7_0806800 


metabolism general 


123.0 


11 


9 


4 


29 


42 


0 


0 


0 


0 


3 


3 


3 


45 


PBANKAJ 22380 


i 


PF3D7_0922200 


metabolism general 


44.8 


25 


16 


0 


4 


20 


1 


2 


3 


0 


0 


0 


3 


23 


PBANKA_082310 


BMC Gen 
biomedc 


PF3D7J 444800 


metabolism general 


40.1 


21 


23 


11 


35 


69 


17 


23 


40 


22 


22 


44 


84 


153 


PBANKAJ 30860 


PF3D7_01 02500 


cell surface & PVM 


181.2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 33270 


comics 2( 
entral.cc 


PF3D7_0202500 


cell surface & PVM 


11.5 


1 


0 


0 


0 


0 


X 


X 


X 


X 


X 


X 


X 


0 


X 


PF3D7_0406200 


cell surface & PVM 


16.7 


3 


9 


34 


17 


60 


X 


X 


X 


X 


X 


X 


X 


60 


X 


)13, 


PF3D7_0423700 


cell surface & PVM 


10.2 


0 


5 


0 


2 


7 


X 


X 


X 


X 


X 


X 


X 


7 


X 


14:2 

471 - 


PF3D7_0707300 


cell surface & PVM 


103.6 


3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_080450 


NJ <-n 
— » Q\ 


PF3D7_0731500 


cell surface & PVM 


174.6 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 33270 




PF3D7_0930300 


cell surface & PVM 


195.7 


77 


0 


3 


0 


3 


0 


0 


0 


0 


0 


0 


0 


3 


PBANKA_083100 


NJ 
Ln 

ON 


PF3D7_1021900 


cell surface & PVM 


267.4 


11 


26 


4 


31 


61 


0 


0 


0 


0 


0 


0 


0 


61 


PBANKAJ 12590/ PBANKAJ 41 230 


PF3D7_1 028700 


cell surface & PVM 


58.1 


0 


0 


1 


11 


12 


0 


0 


0 


0 


0 


0 


0 


12 


PBANKA_051280 




PF3D7_1 033200 


cell surface & PVM 


38.9 


11 


16 


0 


4 


20 


X 


X 


X 


X 


X 


X 


X 


20 


X 




PF3D7J 035300 


cell surface & PVM 


141.1 


0 


0 


1 


0 


1 


X 


X 


X 


X 


X 


X 


X 


1 


X 




PF3D7_1 035700 


cell surface & PVM 


80.3 


10 


0 


0 


0 


0 


X 


X 


X 


X 


X 


X 


X 


0 


X 




PF3D7_1 103500 


cell surface & PVM 


67.3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_094340 




PF3D7_1218800 


cell surface & PVM 


39.6 


0 


0 


2 


5 


7 


0 


0 


0 


0 


0 


0 


0 


7 


PBANKAJ 43440 




PF3D7J 228600 


cell surface & PVM 


86.6 


23 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 44330 




PF3D7_1 335300 


cell surface & PVM 


382.9 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_000220 + 10 more 




PF3D7_0206200 


others 


62.7 


0 


0 


0 


2 


2 


0 


0 


0 


0 


0 


0 


0 


2 


PBANKA_030390 




PF3D7_0207700 


others 


108.7 


11 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_030480 




PF3D7_03 11600 


others 


84.0 


0 


0 


0 


5 


5 


0 


0 


0 


0 


0 


0 


0 


5 


PBANKA_040960 




PF3D7_0320700 


others 


21.0 


0 


0 


0 


8 


8 


0 


0 


0 


0 


0 


0 


0 


8 


PBANKAJ 21 780 




PF3D7_0403800 


others 


83.4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 001 50 




PF3D7_0709900 


others 


274.0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 22050 




PF3D7_0817600 


others 


112.2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_071320 




PF3D7_0930000 


others 


66.8 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_083070 




rY5\J/_ I Uz/oUU 


others 


44.2 


1 7 


9 


14 


1 8 


41 


0 


0 


0 


0 


0 


0 


0 


41 


pdankajq i i yu 




PF3D7_1 033800 


others 


52.3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_051760 




PF3D7_1 103100 


others 


13.0 


2 


2 


5 


4 


11 


1 


0 


1 


0 


0 


0 


1 


12 


PBANKA_094360 




PF3D7_1 207700 


others 


43.4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_060620 


a> 
o 


PF3D7_1 358800 


others 


17.3 


8 


6 


4 


9 


19 


8 


5 


13 


1 


0 


1 


14 


33 


PBANKAJ 13510 


o 

— h 

NJ 



Table 2 Comparison of genes identified by SSH with P. falciparum and P. berghei gametocyte proteome data (Continued) 



PF3D7J 408600 


others 


25.1 


8 


8 


7 


12 


27 


2 


5 


7 


2 


0 


2 


9 


36 


PBANKAJ 03390 


PF3D7J 457000 


others 


47.6 


0 


0 


0 


0 


0 


2 


0 


2 


0 


1 


1 


3 


3 


PBANKAJ 32070 


PF3D7_0202400 


unknown 


141.9 


27 


0 


0 


0 


0 


X 


X 


X 


X 


X 


X 


X 


0 


X 


PF3D7_0411000 


unknown 


179.7 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_061240 


PF3D7_0417400 


unknown 


822.4 


0 


0 


1 


0 


1 


0 


0 


0 


0 


0 


0 


0 


1 


PBANKA_071950 


PF3D7_0508900 


unknown 


370.2 


0 


2 


0 


0 


2 


0 


0 


0 


0 


0 


0 


0 


2 


PBANKAJ 10850 


PF3D7_0601900 


unknown 


14.7 


0 


0 


0 


1 


1 


X 


X 


X 


X 


X 


X 


X 


1 


X 


PF3D7_0704100 


unknown 


425.3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_080180 


PF3D7_0730800.1 


unknown 


32.1 


5 


0 


0 


0 


0 


X 


X 


X 


X 


X 


X 


X 


0 


X 


PF3D7_0904200 


unknown 


34.0 


0 


1 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


1 


PBANKA_041720 


PF3D7J021700 


unknown 


830.2 


0 


0 


0 


1 


1 


1 


0 


1 


0 


0 


0 


1 


2 


PBANKA_050590 


PF3D7_1 023000 


unknown 


48.9 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_050720 


PF3D7_1 024800 


unknown 


171.1 


4 


33 


15 


33 


81 


0 


0 


0 


0 


0 


0 


0 


81 


PBANKA_050900 


PF3D7_1026100 


unknown 


23.9 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_051030 


PF3D7_1 026600 


unknown 


196.5 


0 


0 


0 


0 


0 


X 


X 


X 


X 


X 


X 


X 


0 


X 


PF3D7J 114600 


unknown 


39.6 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKA_093310 


PF3D7_1 124200 


unknown 


49.3 


0 


0 


0 


7 


7 


0 


0 


0 


1 


0 


1 


1 


8 


PBANKA_092410 


PF3D7_1 126700 


unknown 


110.1 


0 


1 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


1 


PBANKA_092170 


PF3D7J 225600 


unknown 


93.1 


0 


1 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


1 


PBANKAJ 44050 


PF3D7_1230100 


unknown 


61.5 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 44480 


PF3D7_1 239400 


unknown 


22.0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 45390 


PF3D7_1316700 


unknown 


73.4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 41 520 


PF3D7J321000 


unknown 


80.6 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 41 930 


PF3D7_1 362600 


unknown 


32.1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 13860 


PF3D7_1 438800 


unknown 


84.0 


0 


0 


0 


2 


2 


0 


0 


0 


0 


0 


0 


0 


2 


PBANKAJ 30270 


PF3D7_1 439600 


unknown 


22.6 


0 


6 


0 


10 


16 


1 


0 


1 


0 


0 


0 


1 


17 


PBANKAJ 30350 


PF3D7J451200 


unknown 


176.2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 31 490 


PF3D7_1 452400 


unknown 


107.5 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


PBANKAJ 31 610 



PF troph, peptide hits for gametocyte-less P. falciparum trophozoites [12,13]; I & II, stage I & II P. berghei gametocyte [12,13]; V, stage V P. falciparum gametocytes [9]; I.II.IV.V, sum of all P. falciparum gametocyte data; 
FG1, peptide hits from highly purified female P. berghei gametocyte sample 1 [22]; FG3, peptide hits from highly purified female P. berghei gametocytes sample 3 [22]; FG(1+3), sum of both; MG1, peptide hits from 
highly purified male P. berghei gametocyte sample 1 [22]; MG3, peptide hits from highly purified male P. berghei gametocyte sample 3 [22]; MG(1+3), sum of both; FG+MG, total number peptide hits from all four P. 
berghei gametocyte samples; total FG.MG.PF, total number of P. falciparum & P. berghei gametocyte peptide hits. 
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Table 3 Genes identified by SSH with P. berghei orthologs showing destabilization in the absence of DOZI and CITH 



P. falciparum 
gene ID 


P 


Protein 


Predicted feature 


MW 


P. berghei ortholog 
gene ID 


DOZI 


CITH 


rrjU/_ I U I I oUU 


7/1 

/4 


l)\ i zz antigen 


KH domain type 1, RNA binding 


1 D 1 .0 


PBANKA_ 


J 21 020 


1 A 

1 .4 


1 1 


rrJU/_U/z4yUU 


U 


kinesin-like protein 




111 A 

zzz.o 


PBANKA_ 


_062240 


1 A 

z.4 


1 .z 


rrjU/_ 1 UoooUU 


U 


plasmepsin VII 


SP, 1 TM, aspartyl protease 


JZ.J 


PBANKA_ 


_05 1 760 


i 1 
z. I 


1 3 


rrjU/_ I 5 I o/UU 


U 


conserved Plasmodium protein 


c D 


73 A 
/D.4 


PBANKA_ 


J 41 520 


1 1 
z. I 


1 A 

1 .4 


PF3D7_1 321 000 


0 


conserved Plasmodium protein 




80.6 


PBANKA_ 


J 41 930 


1.8 


1.6 


PF3D7_0303900 


19 


conserved Plasmodium protein 


SP, phosphatidyl-ethanolamine-binding 
protein PEBP domain 


23.2 


PBANKA_ 


_040250 


2.2 


1.8 


PF3D7_1 103500 


0 


CPW-WPC family protein 


SP, CPW-WPC domain 


67.3 


PBANKA_ 


_094340 


1.3 


2.3 


PF3D7_1 1 15200 


0 


SET domain protein, SET7 


Chromatin regulator 


94.3 


PBANKA_ 


_093250 


2.5 


2.3 


PF3D7J 207700 


0 


blood stage antigen 41-3 precursor 


SP 


43.4 


PBANKA_ 


_060620 


2.5 


2.8 


PF3D7_0904200 


1 


conserved Plasmodium protein 


SP, 1 TM 


34 


PBANKA_ 


_041720 


1.9 


3.2 


PF3D7_1 362600 


0 


conserved Plasmodium protein 


Interacts with splicing factor 


32.1 


PBANKA_ 


J 13860 


3.4 


3.5 



3B subunit 2 



P, total number of peptides identified in P. falciparum and P. berghei gametocytes [9,12,13,22]; DOZI and CITH, log 2 values of transcript abundance WT compared 
to DOZI-KO, or WT compared to CITH KO [20]. 



kinase CLK4 (also termed SRPK1) is expressed in the 
asexual blood stages and gametocytes of P. falciparum 
[45,46] and a gene knock-out of the orthologous protein 
in P. berghei failed to exflagellate upon gametocyte activa- 
tion [47]. The gene set7 encodes for one of the nine SET- 
domain-containing proteins (see above). Transcript of set7 
is destabilized in the absence of DOZI and CITH (compare 
with Table 1), pointing at a translational repression of set7 
in gametocytes and a role in epigenetic control mechanisms 
of the ookinete. 

Pfsl6 is a transmembrane protein of the gametocyte 
PVM [12,27,48,49]. After the PVM has ruptured during 
the egress of the activated gametocyte from the host cell, 
Pfsl6 is not detectable any longer in the sexual stage 
parasites [7] (see below), explaining the decrease of 
transcript following activation. Furthermore, a-tubulin II 
represents a tubulin isoform that is expressed in asexual 
blood stage parasites, in gametocytes and in micro- 
gametes [50-52]. Targeted gene modification studies in 
P. berghei indicated that a-tubulin II plays an import- 
ant role for microgametogenesis [53]. The fact that 
a-tubulin II transcript levels decrease during gameto- 
genesis lets us conclude that the a-tubulin II fila- 
ments needed for the formation of the microgametes 
are already present in the mature gametocytes prior 
to activation. 

The presenilin-like signal peptide peptidase (SPP) was 
previously reported to be involved in merozoite invasion 
and cleavage of the erythrocyte cytoskeletal protein band 
3 [54]. A recent study, however, disagreed with these 
findings, and reported that SPP is an ER-resident pro- 
tease required for growth of the erythrocytic stages [55]. 
SPP is considered a potential drug target of liver and 



blood stage parasites [55-57]. The combined data 
indicate that SPP is present in the liver, blood and 
gametocyte stages, pointing to a role of SPP in multiple 
life-cycle stages of Plasmodium, 

The gene product of PF3D7_08 17600 comprises a 
MAC/PF domain, similar to the MAC/PF domains of 
the previously described plasmodial perforin-like pro- 
teins PPLP1-5 [58] (reviewed in [8]). We therefore 
termed the protein PPLP6 in this study. Furthermore, 
the gene product of PF3D7_1 103500 is a member of a 
family of nine secreted proteins with cysteine-rich CPW- 
WPC domains [59] (reviewed in [60]). The CPW-WPC 
domain is a conserved domain of about 61 residues in 
length that exhibits six well-conserved cysteine residues 
and six well-conserved aromatic sites. The functions and 
life cycle expression patterns of the plasmodial CPW- 
WPC proteins are hitherto not known. 

Three genes were predominantly expressed in the asexual 
blood stages and show a resurgence in transcript expression 
in gametocytes during maturation, i.e. PF3D7_1238900 
(pk2), PF3D7_1 136400 (encoding for a tetratrico peptide 
repeat region; TPR), and PF3D7_1021700. The calcium/cal- 
modulin-dependent protein kinase PK2 was hitherto only 
investigated in the asexual blood stages of P. falciparum 
[61,62] and protein expression in gametocytes has not yet 
been investigated. The gene product of PF3D7_1 136400 
comprises a TPR domain, which is known to mediate 
protein-protein interactions and the assembly of multi- 
protein complexes (reviewed in [63]), but the function of 
the plasmodial TPR domain protein is not yet known. 
While transcript levels of PF3D7_1 136400 increase during 
gametocyte development, a drop following activation was 
observed. 
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(See figure on previous page.) 

Figure 2 Changes in transcript expression levels during the formation and activation of P. falciparum gametocytes. (A) Diagnostic 
RT-PCR demonstrated the presence of stage-specific transcript in the asexual blood stages of the gametocyte-less strain F12 (F12 ABS) as well as 
in immature (stage lll/IV) and mature (stage V) non-activated gametocytes (GC) and in gametocytes at 30 min p.a. in vitro, ama-1, marker for ABS; 
pfccp2, marker for GC. Transcript of gene hdacl was detected in all stages investigated, while no hdac 1 -specific PCR product was detected in 
samples lacking reverse transcriptase (RT). (B) Real-time RT-PCR analysis showed changes of transcript expression of 27 of the SSH-identified 
genes between ABS, immature and mature non-activated GCs and in GCs at 30 min p.a. in vitro. Transcript expression levels were calculated by 
the 2" ACt method; the threshold cycle number (Ct) was normalized with the Ct of the gene encoding for seryl tRNA synthetase as reference. 



For four genes, PF3D7_0818900 {hsp70-l), 
PF3D7_080700 {su a6), PF3D7_0807500 (gapSO), and 
PF3D7_1444800 (encoding for fructose 1,6 bisphosphate 
aldolase, FBPA) expression was high in asexual blood 
stages and decreased during gametocyte development and 
following activation (Figure 2B). SU a6 is a component of 
the plasmodial proteasome, a proteolytic complex com- 
posed of more than 33 SUs that is responsible for the deg- 
radation and recycling of ubiquitinated proteins. As an 
external control we thus investigated the transcript levels 
of another a-ring component, su a5, and revealed similar 
changes in the transcript levels of su aS in asexual blood 
stage parasites and gametocytes before and after activation 
(Figure 2B). 

The genome of P. falciparum encodes for a variety of 
chaperones, including heat shock proteins (HSPs) of the 
HSP70, HSP90 and DnaJ/HSP40 families (reviewed in 
[64]). HSP70-1 was previously described to be located in 
the cytoplasm and nucleus of the parasite blood stages 
[65] as well as in the PV, pointing to a role in protein 
transport to the erythrocyte [66]. Another proposed 
function of HSP70-1 includes the protein trafficking to 
the apicoplast [67,68], indicating that the chaperone has 
several essential functions and is important for multiple 
life-cycle stages. It remains to be elucidated, if other of 
the identified HSPs have more specific functions during 
parasite transmission from human to mosquito. Inte- 
restingly in this context, HSP90 belongs to a family of 
evolutionarily conserved chaperones which has been 
suggested as a capacitor for morphological evolution 
because reduction of its function results in phenotypic 
variation in Drosophila [69]. In addition, HSP90 displays 
multiple roles in stress adaptation and development, 
including spermatogenesis, oogenesis and embryogenesis 
in insects [70]. 

The transmembrane protein GAP50 is part of the IMC 
of the parasite invasive stages, including the ookinete 
[42,71]. Here, it links myosin with the outer membrane 
of the IMC and thus contributes to gliding motility of 
the parasite. Recently, GAP50 was also described as a 
component of the gametocyte IMC [43,72,73], which 
appears to be important for the stability of the crescent- 
shaped cell. Noteworthy, transcript of the P. berghei 
orthologous protein was shown to co-precipitate with 



DOZI and CITH (G.R. Mair, unpublished observations), 
indicating that gapSO is translationally repressed in 
gametocytes. Thus, the release of the translational re- 
pression and onset of protein synthesis during gameto- 
cyte activation might cause the detected decrease in the 
gapSO transcript level. 

FBPA is an enzyme of glycolysis, and in Plasmodium is 
also reported to be part of the motor complex, here 
linking TRAP to actin I [71,74]. Noteworthy in this 
context, during gametocytogenesis, malaria parasites 
shift from glycolysis towards mitochondrial respiration 
[75], which might explain the decrease in FBPA expres- 
sion in gametocytes compared to asexual blood stage 
parasites. 

Transcript expression for five genes, PF3D7_0825800 
{vsp% PF3D7_0215400 (encoding for a WD40 motif), 
PF3D7_1351700 (alv6), PF3D7_1033200 (etrampl0.2), 
and PF3D7_0207700 (sera4), is negligible in game- 
tocytes. The gene PF3D7_08 15800 encodes for a homo- 
log to the yeast vacuolar sorting protein VSP9 with a 
predicted function in protein-protein interactions, while 
PF3D7_0215400 encodes for a protein with a WD40 
motif. WD -repeat proteins are a large family found in all 
eukaryotes, implicated in a variety of functions ranging 
from signal transduction and transcription regulation to 
cell cycle control and apoptosis. The WD40 motifs act 
as a site for protein-protein interaction, and proteins 
containing WD40 repeats are known to serve as plat- 
forms for the assembly of protein complexes or me- 
diators of transient interplay among other proteins 
(reviewed in [76]). WD-repeat proteins of P. falciparum 
have hitherto been described as receptors for protein 
kinase C, and as components of the myosin-driven 
motor complex [77,78]. Furthermore, ALV6 is a member 
of the alveolin family, which comprises seven proteins in 
P. falciparum, associated with the membranous sacs of 
the IMC [79]. The fact that ALV6 expression decreases 
during gametocytogenesis is fairly surprising, conside- 
ring that gametocytes possess an IMC (see above). It has 
to be elucidated, if the expression of any other alveolin 
is up-regulated in gametocytes. 

The early transcribed membrane proteins (ETRAMPs) 
are proteins of the plasmodial PVM [80,81]. In P. falcip- 
arum the proteins were shown to form complexes with 
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the PVM protein exported protein 1 (EXP-1) [82]. To 
date, ETR AMP 10.2 transcripts were found in P. falcip- 
arum trophozoites and the mixed asexual stages and the 
liver stages of P. yoelii [80,81]. Because EXP-1 is also 
present in the gametocyte PVM [7], the expression of 
some of the ETRAMPs in these stages can be expected. 
It thus remains to be elucidated, if the other two SSH- 
identified ETRAMPs, ETRAMP2 and ETRAMP4, might 
play a specific role for gametocytes. Expression of serine 
repeat antigen SERA4 in the asexual blood stages has 
previously been described [83]. The SERA family 
comprises nine proteins with functions in asexual 
blood stage growth and host cell egress (reviewed in 
[8,84]). 

For selected genes we subsequently investigated the 
expression changes on the protein level. Antibodies 
against the proteins PK2, CLK-4, actin I and actin II, 
GAP50, proteasome SU a5, Pfsl6, Pfs25 and PPLP6 
were used to immunolabel the respective proteins in 
samples of P. falciparum F12 asexual blood stages, non- 
activated NF54 gametocytes and gametocytes at 30 min 
p.a. via indirect immunofluorescence assay (Figure 3A). 
Asexual blood stage parasites were highlighted by MSP- 
1 labeling; gametocytes were highlighted by labeling of 
Pfs25 or Pfs230, respectively. The abundance of the 
respective proteins in gametocytes before and after acti- 
vation was quantified by measuring the average fluores- 
cence signal intensity of a total number of 20 plotted 
cells per setting (Figure 3B). A significant up-regulation 
in the expression of PK2, actin II and Pfs25 following 
gametocyte activation was confirmed (Figure 3 A, B). 
While PK2 is also expressed in asexual blood stage para- 
sites, actin II and Pfs25 cannot be detected in these 
stages. On the other hand, CLK4, proteasome SU a5, 
and PPLP6 were detected in asexual blood stage 
parasites and gametocytes during maturation, but the 
proteins were down-regulated in gametocytes at 30 min 
p.a. (Figure 3B). While labeling of CLK4 and SU a5 
revealed a homogenous expression of these proteins in the 
parasite cytoplasm and nucleus, PPLP6-labeling exhibited 
a punctuated expression (Figure 3A). Noteworthy, the 
PPLP6 labeling disappeared after the gametes have fully 
egressed from the enveloping erythrocyte. A role of the 
perforin PPLP1 was assigned to rupturing the PVM 
during host cell egress by the Apicomplexan parasite 
Toxoplasma gondii [85], and a recent study reported the 
involvement of PPLP2 in the lysis of the erythrocyte 
membrane during the exflagellation of male P. berghei 
gametocytes [86]. We thus hypothesize that the plasmo- 
dial PPLP6 might also play a role in the rupture of the 
enveloping membranes of the activated gametocytes, 
explaining its disappearance during gametogenesis. Simi- 
larly, the PVM-based protein Pfsl6 was not present in 
gametocytes at 30 min p.a., since at this time point the 



parasites have egressed from the enveloping erythrocyte 
and the PVM was destroyed [7]. Furthermore, in agree- 
ment with the real-time RT-PCR results, actin I is present 
in both, asexual blood stages and gametocytes, and a 
minor up-regulation in the activated gametocytes was 
detected (Figure 3 A, B). 

Interestingly, we demonstrated a slight up-regulation 
in the protein expression of GAP50 at 30 min p.a., 
although a decrease in transcript levels was described 
during this process. This phenomenon can be explained 
by the fact that translationally repressed gapSO transcript 
is present in the non-activated gametocytes (see above), 
and that the release of repression during activation leads 
to rapid translation and in consequence to a loss in tran- 
script abundance and an increase in protein abundance. 
It is worth mentioning in this context, that we recently 
showed a relocation of GAP 50 from the gametocyte 
IMC to the plasma membrane during gametogenesis. 
Once on the gamete surface, GAP50 binds the human 
complement regulator protein factor H from the blood 
meal, which enables the extracellular parasite to prevent 
lysis by the human complement [43] . 

In summary we demonstrate that the majority of 
genes, which are regulated in their expression levels 
during gametocyte activation, have functions assigned to 
signaling, cell cycle and gene regulation, and proteo- 
stasis, or that they are cytoskeletal or cell surface 
proteins. The regulated expression of proteins involved 
in signaling and cell cycle control is on the one hand 
important for the signaling pathways that result in the 
induction of gametogenesis, once temperature drop and 
XA are perceived by the parasite [47] (reviewed in [6,87]; 
on the other hand gametocyte activation releases the 
parasites from cell cycle arrest, subsequently leading 
to three rounds of rapid DNA replication during micro- 
gametogenesis [88]. 

Furthermore, induction of gametogenesis following 
gametocyte activation results on the one hand in the 
destruction of the gametocyte IMC and relocation of 
cytoskeletal proteins [7,33,34,72,73], on the other hand 
in a drastic turn-over of surface proteins (reviewed in 
[6]). The molecular changes in the composition of cyto- 
skeletal and surface proteins are important for stage 
conversion during gametogenesis and for adaptation of 
the parasite to the insect host and also promote the 
development of the parasite midgut stages, like the 
ookinete. The changes in protein expression during 
parasite transmission from human to insect further 
explain the importance of chaperones and assembly 
proteins. The proteasome on the other hand appears to 
play only a minor role during gametogenesis. We 
hypothesize that proteins which lost their function, once 
the parasite has entered the mosquito vector, are only 
degraded after the gametes have developed and after 
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(See figure on previous page.) 

Figure 3 Changes in protein expression during formation and activation of P. falciparum gametocytes. (A) Immunofluorescence assays, 
using specific antibodies, detected the proteins of interest (in green) in asexual blood stage (ABS) parasites, in gametocytes (GC) and in 
gametocytes at 30 min p.a. (aGC). The parasite stages were highlighted with antibodies against stage-specific markers (in red; MSP1 for ABS; 
Pfs230 or Pfs25 for GCs and aGCs). Nuclei were highlighted by Hoechst nuclear stain (in blue). Bar, 5 urn. (B) Diagram showing the average signal 
intensity of the immunolabeled proteins in gametocytes before and at 30 min p.a. in vitro. Measurements were performed on 20 plotted cells per 
setting via quantitative confocal microscopy. Mean ± SD. *p < 0.05;**p < 0.01; ***p < 0.001 (student's t-test). 



fertilization has occurred, thus after the phase of rapid 
stage conversions has finished. 

Noteworthy, the SSH analysis identified a number of 
unknown proteins with transmembrane domains and/or 
signal peptides, and for six of them, high transcript 
levels were shown in mature and activated gametocytes. 
For three of the genes, expression in the gametocytes 
was higher as compared to asexual blood stage parasites, 
and these genes might encode for proteins important for 
gametocyte development or gametogenesis. 

Conclusions 

The here presented data on the changes in the P. falci- 
parum transcriptome during gametogenesis provide first 
insights into the regulated gene expression occurring 
while the malaria parasite passes the initial phase of 
transmission to the mosquito. Our findings form the 
basis for further studies on genes important for malaria 
transmission and support the identification of proteins 
involved in the development of the mosquito midgut 
stages, which might lead to the discovery of new trans- 
mission blocking targets. 

Methods 

Antibodies 

The following antibodies were used in this study: mouse 
antisera against Pfsl6 and Pfs230 (kindly provided by 
Kim Williamson, Loyola University Chicago), pro- 
teasome SU a5 [35] and actin II [34], as well as rabbit 
antisera against GAP50 (kindly provided by Julian 
Rayner, Sanger Institute England; [89]), Pfs230 (against 
the immunogenic region C) [90], MSP-1 and Pfs25 (BEI 
Resources, Manassas). Mouse antisera against actin I, 
PPLP6, PK2, and CLK4 were generated for this study 
(see below). 

Parasite culture 

Gametocyte-less F12 strain and gametocyte-forming 
NF54 isolate of P. falciparum were cultivated in vitro in 
RPMI 1640 medium complemented with 10% heat- 
inactivated human serum as described [91]. As soon as 
stage II gametocytes started to emerge in the NF54 cul- 
ture, the culture medium was supplemented with 50 
mM N-acetyl glucosamine (GlcNac) for approximately 5 
days to kill the asexual blood stages [92]. The gameto- 
cyte culture was then maintained in normal culture 



medium without GlcNac until immature (stage III and 
IV) or mature stage V gametocytes were harvested and 
enriched by Percoll gradient purification [93]. The ma- 
ture stage V gametocytes were divided into two aliquots, 
and one aliquot was subsequently incubated in game- 
togenesis activating solution (1.67 mg/ml glucose; 8 mg/ml 
NaCl; 1 mg/ml Tris-HCl (pH 8.2) [94], containing 100 uM 
XA for 30 min at RT. 

Construction of a subtracted cDNA library using the SSH 
method 

Total RNA was isolated from Percoll-enriched stage V 
gametocytes before and at 30 min p.a., using Trizol 
reagent (Invitrogen, Karlsruhe, Germany). Subsequently, 
mRNA was isolated from total RNA using the Oligotex 
mRNA Mini Kit (Qiagen, Hilden, Germany) according 
to manufacturers instructions, and the SSH analysis was 
conducted as previously described [21,95] using the PCR- 
Select cDNA Subtraction Kit from Clontech (Mountain 
View, CA, USA), according to the manufacturers pro- 
tocol. Briefly, 220 ng of purified mRNA were reverse- 
transcribed into cDNA using a cDNA synthesis primer; 
subsequently double-stranded cDNA was generated and 
digested with Rsal. The double-stranded cDNA from 
activated gametocytes was ligated in separate aliquots to 
adaptor 1 or adaptor 2R and were denaturated at 98°C for 
90 s and then hybridized at 68°C for 8 h with a 30-fold 
excess of double-stranded cDNA of non-activated game- 
tocytes. Subsequently, both samples were mixed together 
with a 10-fold excess of freshly denatured double-stranded 
cDNA from non-activated gametocytes and hybridized at 
68°C for 16 h. The sample was then subjected to two 
rounds of suppression PCR with PCR-primer 1 and nested 
primers supplied with the kit. PCR amplifications were 
performed in a total volume of 25 \A using a PCR cycler 
(Biometra, Gottingen, Germany) with a heated lid and 
the 5'Prime PCR Extender System (5'Prime, Hamburg, 
Germany). An initial adapter extension at 75°C for 5 min 
was followed by a denaturation step at 93°C for 30 s and 
by 27 cycles of denaturation at 93°C for 15 s, annealing at 
66°C for 30 s, and extension at 72°C for 90 s. A final 
7-min 72°C step was added to allow complete extension 
of the products. The secondary PCR was performed with 
nested primer 1 and 2R on 1 \A of the primary PCR pro- 
ducts for 15 cycles with an initial denaturation step at 
93°C for 1 min, followed by denaturation at 93°C for 15 s, 
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annealing at 68°C for 30 s and extension at 72°C for 90 s, 
plus a final extension step at 72°C for 7 min. Resulting 
PCR products of the secondary subtractive PCR were 
purified using the NucleoSpin Extract II kit (Macherey 
Nagel, Duren, Germany), ligated into the pGEM-T easy 
vector (Promega, Mannheim, Germany) and transformed 
into NEB 5-alpha competent E. coli (New England 
BioLabs, Frankfurt, Germany). The library was plated on 
2xYT agar plates containing 100 (ig/ml ampicillin and 
incubated at 37°C for 16 h. 

Colony PCR, sequencing and computer analysis of cDNA 
sequence data 

A preliminary colony PCR screen on 20 colonies was 
performed with vector-specific primers, i.e. T7-promotor: 
5'- TAATACGACTCACTATAGGG-3' and SP6: 5'-ATTTA 
GGTGACACTATAG- 3' (purchased from Thermo elec- 
tron, Waltham, MA, USA), using the following conditions: 
denaturation at 95°C for 5 min followed by 30 cycles of 
denaturation at 95°C for 15 s, annealing at 43°C for 20 s, 
and extension at 72°C for 60 s. A final 7-min 72°C step 
was added to allow complete extension of the products. 
The colony PCR showed that 75% of clones contained an 
insert in the vector. Subsequently, 288 randomly picked 
clones were screened for the presence of the vector. 
Plasmid isolation of 141 positively screened colonies was 
performed with the Fast Plasmid Mini kit (Eppendorf, 
Hamburg, Germany) and the inserts were sequenced 
(GATC Biotech, Konstanz, Germany; sequences listed in 
the Additional file 1: Table SI). Sequences were used to 
identify similar sequences of the National Center for 
Biotechnology Information databases using BLASTX 
program (BLASTX 2.2.13; http://www.ncbi.nlm.nih.gov/ 
BLAST/) and to predict signal sequences, transmembrane 
regions and features, using the PlasmoDB program (www. 
plasmodb.org) [96]. 

RNA isolation and real-time RT-PCR 

Total RNA was isolated from the mixed asexual F12 
cultures, enriched immature NF54 gametocytes (stage III 
and IV), non- activated mature gametocytes and gameto- 
cytes at 30 min p.a. using the Trizol reagent (Invitrogen) 
according to the manufacturers protocol. RNA prepa- 
rations were treated with RNAse-free DNAse I (Qiagen) 
to remove gDNA contamination, followed by phenol/ 
chloroform extraction and ethanol precipitation. All RNA 
samples had A260/A280 ratios higher than 2.1. Two ug of 
each total RNA sample were used for cDNA synthesis 
using the Superscript III First-Strand Synthesis System 
(Invitrogen), following the manufacturers instructions. 
Controls without reverse transcriptase were used to investi- 
gate potential gDNA contamination by diagnostic PCR 
using hdacl -specific primers. RNA quality was further veri- 
fied for contamination by monitoring transcripts of stage- 



specific genes, i.e. amal and pfccp2, by diagnostic RT-PCR 
(for primer sequences, see Additional file 2: Table S2). 

Primers for real-time RT-PCR were designed using the 
Primer 3 software (http://frodo.wi.mit.edu/primer3/) and 
were initially tested on gDNA in conventional PCR with 
the same conditions subsequently used for real-time 
RT-PCR (for primer sequences, see Additional file 2: 
Table S2). Primers for the reference gene encoding for 
seryl tRNA synthetase were obtained from [38]. The 
PCR products were subsequently separated by agarose 
gel electrophoresis. Primers with high specificity were 
further validated by testing amplification efficiency on 
10-fold dilutions of gDNA using real-time RT-PCR 
followed by melt curve analysis. Primers with low speci- 
ficity, efficiency and poor melt curves were redesigned. 

Real-time RT-PCR measurements were performed 
using the Bio-Rad CFX96 Real-Time Detection System. 
Reactions were prepared in triplicate in a total volume of 
20 ul using the maxima SyBR green qPCR master mix 
(Thermo Scientific, Bonn, Germany) with 20 ng of cDNA 
and primer concentrations of 0.3 uM. The following PCR 
cycling conditions were used: an initial denaturation step 
at 95°C for 10 min, followed by 40 cycles of denaturation 
at 95°C for 15 s and annealing/extension at 60°C for 60 s, 
followed by a final extension at 95°C for 10 s and a final 
melt-curve analysis of 55-95°C. Controls without template 
and without reverse transcriptase were included in all 
real-time RT-PCR experiments. Transcript expression 
levels were calculated by the 2" ACt method [36,37], using 
the endogenous control gene encoding for P. falciparum 
seryl tRNA synthetase (PF3D7_07 17700) as reference 
[38,39]. Control samples lacking reverse transcriptase 
during cDNA synthesis were routinely included dur- 
ing each run and exhibited an average threshold cycle 
number of 39. 



Recombinant protein expression and production of 
mouse antisera 

Recombinant proteins for CLK4, PK2, and actin I were 
expressed as fusion proteins with a GST-tag using the 
pGEX 4 T-l vector (Amersham Biosciences, Freiburg, 
Germany), PPLP6 recombinant proteins were expressed 
as maltose binding protein-tagged fusion proteins using 
vector pIH902 (kindly provided by K. Williamson, 
Chicago; [90]). DNA was amplified by PCR using gene- 
specific primers (for primer sequences, see Additional 
file 2: Table S2). Recombinant proteins were expressed 
in BL21 (DE3) RIL cells according to the manufacturers 
protocol (Invitrogen). Recombinant proteins were isolated 
as either inclusion bodies (PK2, actin I) as described [16] 
or affinity-purified (CLK4 and PPLP6) using glutathione 
sepharose (GE Healthcare, Munich, Germany) and amy- 
lose resin (New England Biolabs), respectively, according 
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to the manusfacturers 'protocols. Specific immune sera 
were generated by immunizing 6-week-old female NMRI 
mice (Charles River Laboratories, Sulzfeld) with 100 \ig 
recombinant protein emulsified in Freunds incomplete 
adjuvant ( Sigma- Aldrich), followed by a boost after 4 
weeks. Sera were collected 10 days after the boost. The 
housing and handling of the animals followed the guide- 
lines of the animal welfare committee of Lower Franconia. 

Indirect immunofluorescence assay and signal strength 
quantification 

Asexual F12 blood stage parasites, non-activated NF54 
gametocytes and gametocytes at 30 min p.a. were air dried 
on glass slides and fixed for 10 min in -80°C methanol. 
For membrane permeabilization and blocking of non- 
specific binding, fixed cells were incubated for 30 min in 
0.01% saponin/0.5% BSA/PBS and 1% neutral goat serum 
( Sigma- Aldrich) in PBS. Preparations were then succes- 
sively incubated for 1.5 h each at 37°C with the primary 
antibody diluted in 0.01% saponin/0.5% BSA/PBS. Binding 
of primary antibody was visualized using fluorescence- 
conjugated goat anti-mouse or anti-rabbit secondary anti- 
bodies (Alexa Fluor 488 or Fluor 596; Molecular Probes, 
Karlsruhe, Germany) diluted in 0.01% saponin/0.5% BSA/ 
PBS. Nuclei were highlighted by incubating the specimens 
with Hoechst nuclear stain 33342 (Molecular Probes) for 
1 min. Cells were mounted with anti-fading solution AF2 
(Citifluor Ltd., Leicester, UK) and sealed with nail varnish. 
Digital images were taken using a Leica TCS SP5 confocal 
laser scanning microscope and processed using Adobe 
Photoshop CS software. For quantitative evaluation of the 
average signal intensity, digital images of 20 randomly 
selected mature and activated gametocytes were taken 
using a Zeiss LSM 510 confocal laser scanning microscope 
under the same optimal laser scanning microscopy 
settings for each antibody. The average fluorescence signal 
intensity of each cell, as indicated by a line edging the 
respective cell, was measured by LSM 510 image software 
and recorded. 
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